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We argue that salient experimental leatures of the magnetic excitations in the SDW phase of iron- 
based superconductors can be understood within an itinerant model. We identify a minimal model 
and use a multi-band RPA treatment of the dynamical spin susceptibility. Weakly-damped spin- 
waves are found near the ordering momentum and it is shown how they dissolve into the particle-hole 
continuum. We show that ellipticity of the electron bands accounts for the anisotropy of the spin 
waves along different crystallographic directions and the gap at the momentum conjugated to the 
ordering one. We argue that our theory agrees well with the neutron scattering data. 

PACS numbers: 74.70.Xa,75.10.Lp,75.30.Ds 

Introduction. The newly discovered iron arsenide superconductors [l| opened a new field in the research on high Tc 
superconductivity. The phase diagram of ferro-pnictide (FP) superconductors is similar to that of layered cuprates 
and contains an antiferromagnetic phase at small dopings and a superconducting phase at larger dopings. However, 
parent compounds of iron-based systems are antiferromagnetic metals rather than Mott insulators, in contrast to the 
cuprates. The electronic structure of parent compounds in the normal state consists of two circular hole pockets of 
unequal size, centered around the F point (0,0), and two elliptic electron pockets centered at (0, ±7r) and (±7r, 0) 
points in the unfolded Brillouin zone (UBZ), which is solely based on the Fe-lattice [l-Q- The dispersions of electron 
and hole bands are significantly nested, i.e. e{j ~ ~£k-i-Q- where Qi is either Qi = (7r,0) or Q2 = (0,7r). Such nesting 
is a boost for antiferromagnetism, and several researchers argued that magnetism is itinerant, and at least partly 
comes from nesting 0-0] • Others argued that magnetism is almost localized and is best described by the Ji — J2 
model Hi. 

Neutron scattering measurements on parent FeAs compounds revealed that magnetic order consists of ferro mag netic 
chains along one crystallographic direction in an Fe plane and antiferromagnetic chains along the other [ol. Il0| i.e., 
the system selects Qi to be either Qi or Q2. Sharp propagating spin- wave like excitations have been observed 
near the ordering momentum (e.g., Qi) up to energies of around 100 meV, with different velocities along the two 
crystallographic directions fs', iTl. At larger energies, excitations become overdamped (ll| . Such stripe order arises 
in the Ji — J2 model of localized spins for J2 < 0.5 Ji, and it was argued that such a model can successfully describe 
some of the experimental data on the spin wave spectra, although, an unusually large in-plane anisotropy of the 
antiferromagnetic exchange between nearest neighbor spins has to be assumed^- However, the same stripe order 
also appears in the itinerant model as a spin-density- wave (SDW) state with either Qi or Q2, and is stabilized by the 
ellipticity of the electron bands and the interactions between the two electron pockets Furthermore, within the 
itinerant model, only one hole and one electron Fermi surface (FS) is involved in the SDW mixing. The other two 
FSs remain intact and give rise to metallic behavior in the SDW phase, even at perfect nesting. 

The issue we address here is whether neutron measurements of magnetic excitations can be explained within the 
itinerant model. This is a crucial test of the itinerant description of magnetism in FPs. A first step in this direction 
was made in Ref. @ , who analyzed spin excitations in the SDW phase of FPs within a two-band model with circular 
hole and electron pockets. We argue that to describe anisotropic magnetic excitations in the SDW state one has to 
consider the model consisting of one circular hole pocket centered around the F-point and two elliptic electron pockets 
centered around the (tt, 0) and (0, vr) points of the unfolded BZ, respectively. Based on this model, we develop a 
multi-band RPA treatment of the dynamical spin susceptibility and compare the results with available experimental 
data. We argue that both the high-energy particle-hole continuum and the low-energy propagating excitations can 
be quantitatively reproduced within our itinerant model. We show that the fact that only one out of two electron 
pockets is involved in the SDW is responsible for the observed anisotropy of the spin wave excitations along different 
crystallographic directions. 

The model. We depart from a 4-band model, use the experimental fact that the two hole FSs around the F-point 
have quite different sizes and assume that one hole band interacts with the electron bands much weaker than the 
other and is thus not important for magnetism. This leads to an effective three band model with one circular hole 
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FIG. 1: (color online) (a) Calculated Fermi surfaces for the three band model. The dashed curves refer to the case of e = 0, 
i.e., complete nesting; Diagrammatic representations of the density-density type interactions (b) and the Dyson equation for 
the RPA spin susceptibility (c) in the three band model. 

pocket at (0,0) (a-band) and two elliptic electron pockets at Qi and Q2 (/3-bands): 

E [^>U^p- + + 4^/5l./32p.] • (1) 

p,CT 

We consider lattice dispersions for all three bands and set Sp = {cos Px + cos py) — /i and e^^ — £0 + 
tp ([1 + e] cos{px + tt) + [1 — e] cos(py)) — ^, — eo + ([1 — e] cos(p:r) + [! + £] cos(pj^ + tt)) — ^. The parameter 
e accounts for the ellipticity of the electron pockets. To make qualitative as well as quantitative comparisons to exper- 
iments, we use the Fermi velocities and the size of the Fermi pockets based on Refs. 0, We obtain ta — 0.85eV, 
tp = — 0.68ey, /i = 1.54ey, eo = O.SleT^, and e = 0.5. For these values, Fermi velocities are Q.beVa for the a-band, 
where a is the Fe — Fe lattice spacing, and Vx = 0.27eVa and Vy = Q.AQeVa along x- and ^/-directions for the /3i-band 
(vice versa for /32, we set ax = ciy — a) . The corresponding Fermi surfaces is shown in Fig. [If a). 

The interacting part of the Hamiltonian contains density-density interactions with small momentum transfer and 
momentum transfers (tt, 0), (0, tt), and (tt, tt). The interactions contributing to SDW formation are (in the terminology 
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of Ref.[i3) 



H4 = Y1 Uial^aP]^^a'PjP2'y'<^Pi'y 
Y [^]p3<T/5lp4<T'"p2'T'apio- + h.c . 



(2) 



For simplicity, we treat both Ui and U3 as constants. Both interactions involve /9i and /32 fermions and in general 
give rise to two SDW OPs Ai oc J^pi^psl^ip+Qn^s-,) with momentum Qi = (7r,0) and A2 oc J^pi'^psl^^p+Q^i'^Sj) 
with momentum Q2 = (0, tt). For perfect nesting and when /3i and /32 fermions do not interact directly, the energy 
depends only on + Ag, i.e., the order is in general a combination of Ai and A2 and the ground state is degenerate. 
However, ellipticity of the electron pockets and a direct interaction between electron bands gives rise to an additional 
aiA^A2 term in the energy, with a positive prefactor ai (Ref. 0)- Then, the energy is minimized when either Ai — 
or A2 = 0, i.e., SDW order is ferromagnetic along one direction and antiferromagnetic along the other, in agreement 
with experiments. Such a state couples a fermions with only one band of /3 fermions, leaving the other band intact, 
i.e., leaving one of the electron FSs unaffected by SDW. 

Without loss of generality we direct Ai along the z-axis and set A2 = 0. The resulting mean-field Hamiltonian can 
be diagonalized by a standard Bogolyubov transformation 



^pa 



/3ip+Qic7 = crivpCpa - Updprr). 



(3) 



where Up,Vp — ^ 1 ± 



'P '^p+Qi 



After the transformation, the Hamiltonian H2 + H4 becomes 



rrdiag _ \^ [ I?c t ^ A-F/H'' H A- /?t 
SDW ~ / . ^p^pa'^va ^ -'^p"pcr"pa' ^ ^p /-'2p(Ti"2p(T 



(4) 



where 



T?c,d 



^p+C+Q.±\/(^p-V+qJ'+4A? 



(5) 



are the new dispersions and the magnitude of the SDW gap is determined self-consistently from Ai = 
-UsDW Ep '^v'^v [/(^p) - /(^p)] > where Usdw = Ui + U3 and f{E) is the Fermi function. 

Magnetic Susceptibility and RPA. We will follow earlier works and use a generalized RPA approach. The 

dynamical susceptibility tensor in the multiband case is defined as 



St. ah/ I •r^\ 
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dre'"'^ 



E 



p,p',7,5,7',5' 

(r.4^(r)W(r)a^,^,(0)6p-_qM'(0))af,,a:;j,, 



(6) 



where is the time ordering operator and quasiparticle operators s, t, a, h are operators from either a or /3i_2 bands. 
Because the unit cell is doubled in the SDW state, the susceptibility is nonzero for q=q' and q=q'+Qi (Ref. 
[H). In terms of Green functions (GF) we have x'*'''''(q, i^^) = -jnT^u Ep <r Gp^(*'^»)Gp+q<x(«w„ + iVL), and 



st.ab / T\\ 



^E.„EpG^t(*'^„)G;,%^(*^„ + ^17), where G;*<,(*w„) = - dr(r.Sp,(r)it^(0))e^-"^ In order 
to compute the RPA susceptibility, we also include all other interactions: the exchange interaction between a and /3 
fermions and the interactions between electron pockets (see Fig. [Hb)). These interactions do not contribute towards 
SDW formation but they do play a role in determining the structure of the spin wave excitations. 

in the RPA approximation are obtained via a Dyson equation 

BPA 



Spin susceptibilities 



Xlr. 



st,ah 
X-lm 



RPA 



st.ab 



st.ab' TTa'h' ,s' t' s't'.ab 



u: 



Xlr. 



(7) 



and the summation over repeated band indices is assumed. The Dyson equation is schematically shown in FigUJc) (de- 
pending on 1,111, the series contain either ladder or bubbles). The solution of Eq.© in matrix form is straightforward 
[xirnlfiPA — Xini(l ~ UimXim)~^ ■ For the single band Hubbard model our results agree with those of Ref. [T5j . 
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FIG. 2: (color online) Calculated imaginary part of the transverse (upper panel) and longitudinal (lower panel) component of 
the physical RPA spin susceptibility as a function of frequency and momentum (plotted along (0, vr) (ji, 0) points of the BZ) 
for e = (a)-(c), and e = 0.5 (b)-(d). The color bars refer to the intensity in units of st./eV. For numerical purposes we set the 
damping constant 5 — l/300eV. 



Results and Comparison to Experiments. We set Usdw ~ 0.52eV for zero ellipticity, to match the experimental 
Tjv ~ 200/^. This yields a T = value of the SDW gap of Ai w 31meV. As the interactions between electron pockets 
are likely smaller than between electron and hole pockets, we set, for definiteness, U5 — U2 — 0.5Ui, Ui = U3, and 
Uji = O.lUj. We verified that the spin wave dispersion around the ordering wave vector does not change in any 
substantial way if we vary these numbers. In particular, when the ellipticity parameter e 7^ 0, we can set Uji = 
and it will not affect the outcome. The results of our calculations are presented in Fig. [2] for T = and for several 
values of the ellipticity parameter e. Consider first the case of circular electron pockets, e = 0, which is for magnetic 
properties the same as full nesting. In this case two out of three FS are completely gapped, if the order parameter 
exceeds the threshold value, which is the case for our A. In Figl2lja) we show the calculated imaginary part of 

the transverse component of the susceptibility, Im [xh (q, fi)]^p^ along the direction (0,7r) to (tt, 0). Because the 

corresponding ground state is degenerate, the excitation spectrum has a Goldstone mode at the ordering momentum 
Qi and another gapless mode at Q2. This result was earlier found in Ref. f^. We also clearly see that the excitations 
near the ordering momentum are propagating up to i7c(Qi) = 2Ai ^ 60meV what is expected because both pockets 
separated by Qi are gapped. Excitations near Q2 are also propagating, but up to smaller energies, rtc{Q2) — Ai, that 
is consistent with the fact that only one of the two pockets separated by Q2 is gapped. At higher energies spin- waves 
enter the continuum and become overdamped Stoner-like excitations, however, still with well defined peaks. The 
difference between nc(Qi) and r2c(Q2) is also seen in Figl2l^c) where we plotted the longitudinal component of the 
susceptibility. As expected, Im[xH (q, f2)]^p^ vanishes below ric(Q). 

At finite e = 0.5, we used a slightly larger Usdw to recover the same Tjy as for e = 0, which in turn leads to a 
slightly larger Ai. The results are presented in panels (b) and (d) of Fig. There are two key effects introduced 
by ellipticity. First, the degeneracy is now broken and the transverse excitations near Q2 acquire a finite gap clearly 
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FIG. 3: (color online) Left panel: Calculated Imx^p^ as a function of frequency and momentum (plotted along (0, 27r) — 
(0, tt) — )■ (tt, tt) directions of the BZ) for e = 0.5. The crosses and circles are the measured points taken from Refs. [llj and 
[31, respectively. Right panel: Calculated constant energy cuts of Imx"' (q, f2) along (0,0) (tt, 0) — > (27r,0) direction (solid 
curves) for 25 (a), 48 (b), 64 (c), 100 (d), 115 (e), 135 (f), and 144 (g) (all in meV). The experimental data are taken from 
Ref.@]. The intensity maxima have been adjusted to the experimental values. 



visible in Figl2fc). Second, spin-wave excitations near Qi have a finite Landau damping down to il = because 
the SDW order no longer completely gaps the hole and electron FSs separated by Qi. At the same time, spin- wave 
excitations are still clearly visible at low energies and do not become overdamped. The reasoning is that SDW 
coherence factors suppress the scattering between a and /3i fermions and reduce the Landau damping around the 
ordering momentum Qi from fi/jq — Qi| to r2|q — Qi| what for linear dispersion is of the same order as and 
(q— Qi)^. Within the one-band model, this effect has been discussed in Once exceeds r^c(Q2), the coherence 
factors no longer screen Landau damping, and spin-waves become Stoner-like, overdamped excitations. The scattering 
between c and /?2 fermions is not suppressed by the SDW coherence factors because Q2 is not the ordering momentum, 
therefore, Landau damping makes the transverse excitations near Q2 overdamped immediately above the gap. Note, 
that there is also a tendency towards incommensuration near Q2. Longitudinal excitations are still gapped up to 
flciQ) and become Stoner-like at higher energies, as is clearly seen in Fig. [^d). 

Despite the fact that our model is indeed a simplification of the actual 5-band model for FPs, it reproduces quite 
well the experimental data for the magnetic excitations. First, the excitations near the "wrong" momentum Q2 are 
gapped up to about 50meV as seen in Fig[2I^c). Second, measured excitations near Qi are propagating spin-waves 
up to about 100 meV, with a finite but small damping. In the left panel of Fig[3]we compare our dispersion with the 
experimental data [ll| . We emphasize that both the measured and the calculated spectra have different velocities 
along different crystallographic directions. In our theory, it is a consequence of the non-zero ellipticity and the fact 
that only one electron pocket is involved in the SDW formation. Remarkably, for e = 0.5 the anisotropy is the same as 
in the experimental data. Another verifiable theory prediction for this range is that the width of the spin-wave peak 
should scale as fl^. Third, we find that excitations are still visible even when the spin- waves enter the continuum. 
In the right panel of Figl3] we compare the calculated and the measured Imx^ (9, ^) along (0, 0) — » (tt, 0) (27r, 0) 
direction for different Q.. We see that the agreement is quite good for all frequencies even above ilc [17| . 

Conclusions In this paper we analyzed the structure of the magnetic excitations in the magnetically ordered state of 
Fe-pnictides. We used a multi-band itinerant model and developed a multi-band RPA treatment of the longitudinal and 
transverse components of the dynamical spin susceptibility. We found weakly-damped spin-waves near the ordering 
momentum and showed how they dissolve into the particle-hole continuum above an energy which scales with the 
SDW order parameter. For perfect nesting between electron and hole bands the SDW state has an extra degeneracy, 
and we found an extra gapless mode at the momentum different from the ordering one. When ellipticity of the electron 
bands is included, the degeneracy is lifted and the extra mode becomes gapped. Ellipticity of the electron pockets 
also naturally accounts for the anisotropy of the spin waves along different crystallographic directions. We argue that 
our theory agrees well with the available neutron scattering data. 
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